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Abstract 

Using a family of modified WeibuU distributions, encompassing both sub-exponentials and super-exponentials, 
to parameterize the marginal distributions of asset returns and their natural multivariate generalizations, 
we give exact formulas for the tails and for the moments and cumulants of the distribution of returns 
of a portfolio make of arbitrary compositions of these assets. Using combinatorial and hypergeometric 
functions, we are in particular able to extend previous results to the case where the exponents of the 
Weibull distributions are different from asset to asset and in the presence of dependence between assets. 
We treat in details the problem of risk minimization using two different measures of risks (cumulants 
and value-at-risk) for a portfolio made of two assets and compare the theoretical predictions with direct 
empirical data. While good agreement is found, the remaining discrepancy between theory and data 
stems from the deviations from the Weibull parameterization for small returns. Our extended formulas 
enable us to determine analytically the conditions under which it is possible to "have your cake and eat 
it too", i.e., to construct a portfolio with both larger return and smaller "large risks". 



1 Introduction 

The determination of the risks and returns associated with a given portfolio constituted of assets is com- 
pletely embedded in the knowledge of their multivariate distribution of returns. Indeed, the dependence 
between random variables is completely described by their joint distribution. This remark entails the two 
major problems of portfolio theory: 1) determine the multivariate distribution function of asset returns; 2) 
derive from it useful measures of portfolio risks and use them to analyze and optimize portfolios. Here, 
we address them both by extending the new approach of [^] in terms of a class of subexponential and 
superexponential multivariate distributions. 

In the standard Gaussian framework, the multivariate distribution takes the form of an exponential of minus 
a quadratic form X'Q^^X, where X is the unicolumn of asset returns and Q, is their covariance matrix. 
The beauty and simplicity of the Gaussian case is that the essentially impossible task of determining a 
large multidimensional function is collapsed into the very much simpler one of calculating the N{N + 
l)/2 elements of the symmetric covariance matrix. Risk is then uniquely and completely embodied by the 
variance of the portfolio return, which is easily determined from the covariance matrix. This is the basis of 
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Markovitz's portfolio theory [17] and of the CAPM (see for instance [18|). In this framework, increasing 



return comes with risk and return may remunerate a larger risk. 

However, the variance (volatility) of portfolio returns provides at best a limited quantification of incurred 
risks, as the empirical distributions of returns have "fat tails" |8|] and the dependences between assets are 
only imperfectly accounted for by the co variance matrix [|l^]. Value-at-Risk [[IT]] and other measures of risks 
^ ^] have been developed to account for the larger moves allowed by non-Gaussian distributions 
and nonlinear correlations. 

In section 2, we present our parameterization of the multivariate distribution of returns based on two steps: 

(i) the projection of the empirical marginal distributions onto Gaussian laws via nonlinear mappings; (ii) 
the use of an entropy maximization to construct the corresponding most parsimonious representation of the 
multivariate distribution. We show in particular that this construction amounts to use Gaussian copulas. 
Empirical tests are given which show that this assumption is a very good approxition. 

Section 3 offers a specific parameterization of marginal distributions in terms of so-called modified Weibull 
distributions, which are essentially exponential of minus a power law. This family of distribution contains 
both sub-exponential and super-exponentials, inclusing the Gaussian law as a special case. Notwithstanding 
their possible fat-tail nature, all their moments and cumulants are finite and can be calculated. We present 
empirical calibration of the two key parameters of the modified Weibull distribution, namely the exponent c 
and the characteristic scale x- 

Section 4 uses the multivariate construction based on (i) the modified Weibull marginal distributions and 

(ii) the Gaussian copula to derive the asymptotic analytical form of the tail of the distribution of returns 
of a portfolio composed of an arbitrary combination of these assets. We show that, in the case where 
individual asset returns have the asymptotic tail ex]3[—{xi/xiy] with the same exponent c, then the tail of 
the distribution of portfolio return S is of the same form ex.p[—{S/xy], with the same exponent c but with 
a characteristic scale x taking different functional forms depending on the value of c and the strength of the 
dependence between the assets, x is also of course a function of the asset weights in the portfolio. These 
results allow one to estimate the value-at-risk (VaR) in this non-Gaussian nonlinear dependence framework. 

Section 5 provides the analytical expressions of the cumulants of the distribution of portfolio returns. Cumu- 
lants are of interest because they are natural measures of risks, the higher their order, the more weight being 
given to large risks. Recall that the second cumulant is nothing but the variance, the normalized third-order 
(resp. fourth-order) cumulant is the skewness (resp. excess kurtosis). This section provides the most general 
formulas for any possible positive values of the exponents c of the asset return distiibutions, generalizing 



generously previous results [26]. In the case of dependent assets, we give the rather cumbersome explicit 
formulas only for the case of portfolios of two assets. Similar more cumbersome expressions hold for the 
general case of N assets. This section also offers empirical tests comparing the direct numerical evalua- 
tion of the cumulants of financial time series to the values predicted from our analytical formulas using the 
exponents c and characteristic scales x calibrated previously on the same data. Good consistency is found. 

Section 6 uses these two sets of results to offer a first approach to portfolio optimization. We use two ap- 
proaches, one based on the asymptotic form of the tail of the distribution of portfolio returns, the other based 
on the expression of the cumulants of the distribution of the portfolio returns. In both cases, we show how to 
generalize the concept of an efficient frontier, initially introduced in the mean-variance space. We extend it 
using the different measures of risks captured by the tail of the distribution and by the cumulants. The main 
novel result is an analytical understanding of the conditions under which it is possible to simultaneously in- 
crease the portfolio return and decreases its large risks quantified by large-order cumulants. It thus appears 
that the multidimensional nature of risks allows one to break the stalemate of no better return without more 
risks. Section 7 concludes. 
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Before proceeding with the presentation of our results, we set the notations to derive the basic problem 
addressed in this paper, namely to study the distribution of the sum of weighted random variables with 
arbitrary marginal distributions and dependence. Consider a portfolio with shares of asset i of price Pi{0) 
at time t = whose initial wealth is 



N 



w{o) = Y,nm{o) ■ (1) 



i=l 

A time r later, the wealth has become W{t) = J2iLi '^iVii'T') and the wealth variation is 

5rW = W{t) - W{Q) = ^nmi0 f'^^' ' = W{0) ^wMt,T), (2) 

i=i P'^^' i=i 

where 

niPi{0) 

Wi = — (3) 

is the fraction in capital invested in the ith asset at time and the return rj(t, r) between time t — t and t of 
asset i is defined as: 

r.(t.r)^ "■">-"•"--' . (4) 

Pi{t-T) 

Using the definition (^), this justifies us to write the return Sr of the portfolio over a time interval r as the 
weighted sum of the returns rj(r) of the assets i = 1, ...,N over the time interval r 

= WW) ^ h"^' "'^^^ ' ^ ^ 

In the sequel, we shall thus consider asset returns as the fundamental variables (denoted Xi or Xi in the 
sequel) and study their aggregation properties, namely how the distribution of portfolio return equal to their 
weighted sum derives for their multivariable distribution. We shall consider a single time scale r which 
can be chosen arbitrarily, say equal to one day. We shall thus drop the dependence on r, understanding 
implicitely that all our results hold for returns estimated over time step r. 



2 Estimation of the joint probability distribution of returns of several assets 
2.1 A brief exposition and justification of the method 



We will use the method of determination of multivariate distributions introduced by Karlen [12] and Somette 



et al. [|26|]. This method consists in two steps: (i) transform each return x into a Gaussian variable ?/ by a 
nonlinear monotonous increasing mapping; (ii) use the principle of entropy maximization to construct the 
corresponding multivariate distribution of the transformed variables y. 

The first concern to address before going any further is whether the nonlinear transformation, which is in 
principle different for each asset return, conserves the structure of the dependence. In what sense is the 
dependence between the transformed variables y the same as the dependence between the asset returns x? It 
turns out that the notion of "copulas" provides a general and rigorous answer which justifies our procedure 

For completeness, we briefly recall the definition of a copula (for further details about the concept of copula 



see [20]). A function C : [0, 1]" — > [0, 1] is a n-copula if it enjoys the following properties : 
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• Vug [0,1], •••,!) = u, 

• Vuj G [0, 1], C(ni, • • • , Un) = if at least one of the Uj equals zero , 

• C is grounded and n-increasing i.e the C-volume of every boxes whose vertices lie in [0, 1]" is posi- 
tive. 

Skar's Theorem then states that, given an n-dimensional distribution function F with continuous marginal 
distributions Fi, • • • , F„, there exists a unique n-copula C : [0, 1]" — > [0, 1] such that : 

F{xi, ■ ■ ■ ,Xn) = C{Fi{xi), ■ ■ ■ , FniXn)) ■ (6) 

This elegant result shows that the study of the dependence of random variables can be performed inde- 
pendently of the behavior of the marginal distributions. Moreover, the following result shows that copulas 
are intrinsic measures of dependence. Consider n continous random variables Xi, - ■ ■ , X„ with copula C. 
Then, if gi{Xi), ■ ■ ■ ,gn{Xn) are strictly increasing on the ranges of Xi, - ■ ■ ,Xn, the random variables 



Yi = gi{Xi), ■ ■ ■ ,Yn = gn{Xn) havc exactly the same copula C [14]. The copula is thus invariant under 
strictly increasing tranformation of the variables. This provides a powerful way of studying scale-invariant 
measures of associations. It is also a natural starting point for construction of multivariate distributions and 
provides the theorical justification of the method of determination of mutivariate distributions that we will 
use in the sequel. 



2.2 Transformation of an arbitrary random variable into a Gaussian variable 

Let us consider the return X, taken as a random variable characterized by the probability density p{x). The 
transformation y{x) which obtains a standard normal variable y from x is determined by the conservation 
of probability: 

1 _^ 

p{x)dx = e 2 dy . (1) 



'2ir 

Integrating this equation from —oo and x, we obtain: 

1 



F(x] 



(8) 



1 +erf 

where F{x) is the cumulative distribution of X: 

F{x) = / dx'p{x') . (9) 



This leads to the following transformation y{xy. 

y = V2erf"^(2F(x)-l), (10) 

which is obvously an increasing function of X as required for the application of the invariance property of 
the copula stated in the previous section. An illustration of the nonlinear transformation ( p^ is shown in 
figure [!} Note that it does not require any special hypothesis on the probability density X, apart from being 
non-degenerate. 



In the case where the pdf of X has only one maximum, we may use a simpler expression equivalent to (|10|). 
Such a pdf can be written under the so-called Von Mises parametrization [|5|] : 

= C^P=e-i/(-) , (11) 
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where C is a constant of normalization. For —>■ when |x| —>■ +00, the pdf has a "fat tail," i.e., it 

decays slower than a Gaussian at large 

Let us now define the change of variable 

y = sgn{x)y/\f{x)\ . (12) 



Using the relationship p{y) = p{x)^, we get: 



1 -2 



P{y) = -i^e 2 . (13) 



It is important to stress the presence of the sign function sgn{x) in equation (|12[), which is essential in order 
to correctly quantify dependences between random variables. This transformation (|l^) is equivalent to ((K 
but of a simpler implementation and will be used in the sequel. 



2.3 Determination of the joint distribution : maximum entropy and Gaussian copula 

Let us now consider N random variables Xi with marginal distributions Pi{xi). Using the transformation 
([lO|), we define N standard normal variables Yi. If these variables were independent, their joint distribu- 
tion would simply be the product of the marginal distributions. In many situations, the variables are not 
independent and it is necessary to study their dependence. 

The simplest approach is to construct their covariance matrix. Applied to the variables Yi, we are certain that 
the covariance matrix exists and is well-defined since their marginal distributions are Gaussian. In contrast, 
this is not ensured for the variables Xi. Indeed, in many situations in nature, in economy, finance and in 
social sciences, pdf 's are found to have power law tails ~ ^^r^ for large If /i < 2, the variance and the 
covariances can not be defined. If 2 < /i < 4, the variance and the covariances exit in principle but their 
sample estimators converge poorly. 

We thus define the covariance matrix: 

V = E[yy^] , (14) 

where y is the vector of variables Yi and the operator E[-] represents the mathematical expectation. A clas- 
sical result of information theory [^] tells us that, given the covariance matrix V, the best joint distribution 
(in the sense of entropy maximization) of the N variables Yi is the multivariate Gaussian: 

P(y) = ,^ , exp ( --yV'V 1 ■ (15) 

Indeed, this distribution implies the minimum additional information or assumption, given the covariance 
matrix. 

Using the joint distribution of the variables Yi, we obtain the joint distribution of the variables Xi. 

dyi 



P(x) = P(y) 

Since 



dxj 



(16) 



dy: 

we get 



dxj 



1 0,2 



dyi 



dxj 



/27Tpj{xj)e-2y^6ij , (17) 

N 

(2^)^/2 J] p,(x,)ei^? . (18) 



i=l 
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This finally yields 

1 ( \ \ ^ 

As expected, if the variables are independent, V = I, and P(x) becomes the product of the marginal 
distributions of the variables Xi. 

Let F(x) denote the cumulative distribution function of x and Fi{xi) their marginal distributions. The 
copula C is then such that 

F{xi,---,Xn)=C{Fi{xi),---,Fn{Xn)) • (20) 

Differentiating with respect to xi , • • • , a:Ar leads to 

P(xi, • • • , = ^-^(^i'--^'^") = c(Fi(:ei), • • • , F„(x„)) TTp^C^O , (21) 

OXi ■ ■ ■ OXn 

1=1 

where 

9C(ui,---,MAr) 

ci^i, • • • ,nAr) = ^ (22) 

OUi ■ ■ ■ OUn 

is the density of the copula C. 

Comparing (22) with ([l9|), the density of the copula is given in the present case by 

c(ni,---,ujv) = -^==exp (^-^y{^)(y~^ -/)y(u) ^ , (23) 

which is the "Gaussian copula" with covariance matrix V. This result clarifies and justifies the method of 
[ ^ ] by showing that it essentially amounts to assume arbitrary marginal distributions with Gaussian copulas. 
Note that the Gaussian copula results directly from the choice of maximizing the Shannon entropy. This is 
not unexpected in analogy with the standard result that the Gaussian law is minimizing the Shannon entropy 
at fixed given variance. If we were to extend this formulation by considering more general expressions of 
the entropy, such that Tsallis entropy [27], we would have found other copulas. 



2.4 Empirical test of the Gaussian copula assumption 

We now present preliminary tests of the hypothesis of Gaussian copulas between returns of financial assets. 
Testing the gaussian copula hypothesis is a delicate task. A priori, two standard methods can be proposed. 

The first one consists in using that Gaussian variables are stable in distribution under addition. Thus, a 
(quantile-quantile or Q — Q) plot of the cumulative distribution of the sum yi + ■ ■ ■ + Vp versus the cumu- 
lative Normal distribution with the same estimated variance should give a straight line in order to qualify 
a multivariate Gaussian distribution (for the transformed y variables). Such tests on empirical data are 
presented in figures 

The second test amounts to estimating the covariance matrix V of the sample we consider. This step is 
simple since, for fast decaying pdf 's, robust estimators of the covariance matrix are available. We can then 
estimate the distribution of the variable = y*V^-'^y. It is well known that follows a distribution 
if y is a Gaussian random vector. Again, the empirical cumulative distribution of versus the cumulative 
distribution of should give a straight line in order to qualify a multivariate Gaussian distribution (for the 
transformed y variables). Such tests on empirical data are presented in figures 
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First, one can observe that the Gaussian copula hypothesis appears better for stocks than for currencies. Note 
also that the test of aggregation seems systematically more in favor of the Gaussian copula hypothesis than 
is the test, maybe due to its smaller sensitivity. Nonetheless, the very good performance of the Gaussian 
hypothesis under the aggregation test bears good news for a porfolio theory based on it, since by definition 
a portfolio corresponds to asset aggregation. Even if sums of the transformed returns are not equivalent to 
sums of returns (as we shall see in the sequel), such sums qualify the collective behavior whose properties 
are controlled by the copula. 

Notwithstanding some deviations from linearity in figures it appears that, for our purpose of developing 
a generalized portfolio theory, the Gaussian copula hypothesis is a good approximation. A more systematic 
test of this goodness of fit requires the quantification of a confidence level, for instance using the Kolmogorov 
test, that would allow us to accept or reject the Gaussian copula hypothesis. However, we have encountered 
a practical problem to implement this test, as the same sample is used both for the empirical calibration of 
the covariance and for the test itself. A bootstrap method is thus necessary in order to obtain a real measure 
of the departure from Gaussianity and thus to decide wether or not these deviations are significant. 



3 Choice of an exponential family to parameterize the marginal distribu- 
tions 

3.1 The modified WeibuU distributions 

We now apply these constructions to a class of distributions with fat tails, that have been found to provide 
a convenient and flexible parameterization of many phenomena found in nature and in the social sciences 
[p^]. These so-called stretched exponential distributions can be seen to be general forms of the extreme tails 
of product of random variables [^. 



Following [g6p, we postulate the following marginal probability distributions of returns: 



1 c -(\iiY 

p{x) = -^ — \x\2-^e , (24) 

where c and x are the two key parameters. A more general parameterization taking into account a possible 
asymmetry between negative and positive returns (thus leading to possible non-zero average return) is 

pix) = ^^-^|x|^-^e~(^) ifx>0 (25) 

p{x) = ^^-^|x|^"^e"(^) ifx<0. (26) 

Thes expressions are close to the WeibuU distribution, with the addition of a power law prefactor to the 
exponential such that the Gaussian law is retrieved for c = 2. Following [26, |5|, |2|], we call ( ^4| ) the 



modified WeibuU distribution. For c < 1, the pdf is a stretched exponential, also called sub-exponential. 
The exponent c determines the shape of the distribution, fatter than an exponential if c < 1. The parameter 
X controls the scale or characteristic width of the distribution. It plays a role analogous to the standard 



deviation of the Gaussian law. See chapter 6 of [24] for a recent review on maximum likelihood and other 



estimators of such generalized Weibull distributions. 



7 



3.2 Transformation of the modified Weibull pdf into a Gaussian Law 



One advantage of the class of distributions (24) is that the transformation into a Gaussian is particularly 
simple. Indeed, the expression (|24|) is of the form (ITT]) with 



X 



\x\ 



Applying the change of variable (O) which reads 



(27) 



yi = sgn{x. 



Xi 



leads automatically to a Gaussian distribution. 

These variables Yi then allow us to obtain the covariance matrix V : 

^ - \ 2 / X 



^ij = tY^ sgn{xiXj) ( — ^ ( 



n=l 



Xi 



Xj 



and thus the multivariate distributions -P(y) and P(x 

1 



P{xi, ■ ■ ■ ,Xn) 



2N^N/2^±L c/2 

1=1 Ai 



n 



■ exp 



c/2 



c/2 



(28) 



(29) 



(30) 



Similar transforms hold, mutatis mutandis, for the asymmetric case. 



3.3 Empirical tests and estimated parameters 



In order to test the validity of our assumption, we have studied a large basket of financial assets including 
currencies and stocks. As an example, we present in figures || to |T2| typical log-log plot of the transformed 
return variable Y versus the return variable X for a certain number of assets. If our assumption was right, 
we should observe a single straight line whose slope is given by c/2. In contrast, we observe in general 
two approximately linear regimes separated by a cross-over. This means that the marginal distribution of 
returns can be approximated by two modified Weibull distributions, one for small returns which is close to a 
Gaussian law and one for large returns with a fat tail. Each regime is depicted by its corresponding straight 
line in the graphs. The exponents c and the scale factors x for the different assets we have studied are given 
in tables |l]for currencies and|2|for stocks. The coefficients within brackets are the coefficients estimated for 
small returns while the non-bracketed coefficients correspond to the second fat tail regime. 

The first point to note is the difference between currencies and stocks. For small as well as for large returns, 
the exponents c_ and c+ for currencies (excepted Poland and Thailand) are all close to each other. Additional 
tests are required to establish whether their relatively small differences are statistically significant. Similarly, 
the scale factors are also comparable. In contrast, many stocks exhibit a large asymmetric behavior for large 
returns with c+ — c_ > 0.5 in about one-half of the investigated stocks. This means that the tails of the large 
negative returns ("crashes") are often much fatter than those of the large positive returns ("rallies"). 

The second important point is that, for small returns, many stocks have an exponent (c+) ^ (c_) ~ 2 and 
thus have a behavior not far from a pure Gaussian, while the average exponent for currencies is about 1.5 
in the same "small return" regime. Therefore, even for small returns, currencies exhibit a strong departure 
from Gaussian behavior. 
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In conclusion, this empirical study shows that the modified WeibuU parameterization, although not exact on 
the entire range of variation of the returns X, remains consistent within each of the two regimes of small 
versus large returns, with a sharp transition between them. It seems especially relevant in the tails of the 
return distributions, on which we shall focus our attention next. 

4 Asymptotic estimation of the Value-at-Risk 

Here, we examine theoretically the tail of the distribution of returns of portfolios constituted of assets with 
distributions characterized by the modified Weibull distributions. Two distinct situations can occur. 

• The tail exponents c of the distributions of asset returns are different from asset to asset. This is the 
most general case and at the same time the simplest. When the assets have different exponent c, the 
asymptotic tails of the portfolio return distribution are dominated by the asset with the heaviest tails. 
The largest risks of the portfolio are thus controlled by the single most risky asset characterized by 
the smallest exponent c. Such extreme risk cannot be diversified away. The best strategy focused on 
minimizing the extreme risks in such a case consists in holding only the asset with the thinnest tail, 
i.e., with the large exponent c. 

• All assets in the portfolio have the same tail exponent c. This case is the most interesting and chal- 
lenging as we shall see. We now present general asymptotic results for this case. 

4.1 Portfolio made of assets with the same exponent c > 1 
4. 1. 1 Case of independent assets 

Consider N assets caracterized by their returns Xi, i = {1, 2, • • • , N}, with joint distribution P(x). We 
have seen in the first section that the portfolio return S can be written as 

N 



where the weights Wi are real coefficients. The variables {X^} may be subjected to several constraints. 
Then, the probability density Ps of the random variable S is given by 



where 5 is the Dirac function. Our purpose here is to evaluate its asymptotic behavior in the case where 




(31) 



i=l 




(32) 



P{xi,--- 



xn) = Pi{xi) • • • Pn{xn) ■ 



(33) 



Using a saddle point approximation (see appendix ^), we show that for large S: 




e 



( 



X 



s 



) 



c 



(34) 



where x is given by 



N 




(35) 



i=l 
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Thus, for a large enough VaR or equivalently for a small enough loss probability, 

N _ _ .-VaR 



Fr{S < -VaR} 



N+l 
C 2 



N+l 



x2—{c-l) 



I 2 ^2{c-l) 



i=l 



,5 I MV-%'(^)\36) 



N-1 
C 2 



N 



N+l 



2— (c- 1)" 



c-1 I 1 



2 "^2{c-l) p 



X 



c-1 /Fai? 



(37) 



where r(-, •) is the incomplet Gamma function. Using an asymptotic expansion of the incomplet gamma 
function, this finally yields 



Pr{5 < -VaR} 



JV-l 
C 2 



N 



N+l 



2— (c- 1)" 



c-1 I 1 



2 ' 2{c-l) 



i=l 



exp 


[-(^) 


c- 







(38) 



Our result here recover the previously announced [ ]26[ ] asymptotic form e \ ^ for the tail of the distribu- 
tion but corrects an error in the calculation of the characteristic scale x- 



4.1.2 Case of dependent assets 



We now consider the case of dependent assets with pdf given by equation ( |30| ) or its asymmetric expres- 
sion. The calculation performed in appendix ^ leads to the following asymptotic expression for the density 
function P{S) for large portfolio returns S : 



P{S) 0^ ( ' e~(C 



and the following asymptotic VaR : 



pr{5<-y„ii).c"r(£^.(^)'). 

where C** is a term independent of S, whose expression is given appendix ^ and x is such that 



-1 c/2 c/2 ' 



and the {ctj} are solution of 



c/2 c/2 Y^TZ-l c/2 c/2 



(39) 



(40) 



(41) 



(42) 



Obviously, if the assets are independent, V^j ^ = bij, the cjj's are solution of 



_c-l 



AwiXi, Vz G {1, • • • , A^} 



(43) 
(44) 



X is then simply given by equation (| 



10 



4.2 Portfolio made of assets with the same exponent c < 1 



In the case of independent assets whose tails of distributions are fatter than an exponential, the tail of the 
distribution of price variations of the portfoUo is given by 



P{S) ~ exp 

with 



1^ 

X 



(45) 



X = max {wiXir ■ ■ ,WnXn} ■ (46) 
ie{i,..N} 

This single variable x controls the asymptotic behavior of the tail of the distribution. 



Using the equation (|45|), the Value-at-Risk is solution of 



-VaR 

Pr{5 < -VaR} ~ / dS e \ ^ J , (47) 



oo 



I.e : 

Ft{S < -VaR} r^r (^^,YaR''^ . (48) 

In the case of correlated assets, we have not yet found an asymptotic formula. But, using the standard result 
that a sum of stretched exponentially distributed random variables behaves, in the regime of large deviation, 
like the variable whose deviation is the most extreme, we expect the previous result to remain true. 

4.3 Summary 

In summary, the density function P{S) behaves for large S like 

P{S) ~ e"(^) , (49) 

where 



X 



^^{wiXi)"^'^ ) if c > 1 and independent assets, (50) 



X = r if c > 1 and dependent assets, (51) 

/ ^ c/2 c/2\~ 

( \^''^ 

m.ax {w\^^,--- ,w\^^}\ ifc<l. (52) 
Vie{i,..Af} ) 



X 



Note that the expression ( |50| ) for c > 1 retrieves the result ( |52| ) for c < 1 by taking the limit c ^ 1+, 



showing the continuity of the formulas. 

Let us translate these formulas in intuitive form. For this, we define a value-at-risk (VaR) VaR* which is 
such that its typical frequency is I/Tq. Tq is by definition the typical recurrence time of a loss larger than 
VaR* . In our present example, we take Tq equals 1 year for example, i.e., VaR* is the typical annual shock 
or crash. The expression (p9|) then allows us to predict the recurrence time T of a VaR equal to a times this 
reference value VaR": 

T\ , „ . / VaR* ^ " 



ln( — ) ~ (a^- 1) ( — — ) +C)(lna). (53) 
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Figure ^ shows In ^ versus a. For fixed a, T increases more and more slowly when the exponent c 
decreases. This quantifies our expectation that large losses occur more frequently for sub-exponential dis- 
tributions than for super-exponential ones. 



5 Cumulant expansion of the portfolio return distribution 
5.1 link between moments and cumulants 

Before deriving the main result of this section, we recall a standard relation between moments and cumulants 
that we need below. 

The moments M„ of the distribution P are defined by 

where P is the characteristic function, i.e., the Fourier transform of P : 

/ + 00 
dS P(S)e^^^ . (55) 
-oo 

Similarly, the cumulants Cn are given by 

P{k) = exp . (56) 



\n=l 



Differentiating n times the equation 



nl / — ' nl 

\n=0 / n=l 



we obtain the following recurrence relations between the moments and the cumulants : 



n—l / ^ 

n—\ 



M„ = ^ ]MpCn-p, (58) 



p=0 ^ ^ 

Cn = Mn-Y[^_^jCpMn-p . (59) 
In the sequel, we will first evaluate the moments, which turns out to be easier, and then using eq (59) we 



will be able to calculate the cumulants. Cumulants are indeed the natural objects to quantify risks: as 
seen from their definition (|5^, cumulants of order larger than 2 quantify deviation from the Gaussian law, 
and thus large risks beyond the variance (equal to the second-order cumulant). More importantly, they are 
invariant with respect to translations or change of a return of reference and are thus appropriate measures of 
fluctuations. 
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5.2 Symmetric assets 



We start with the expression (32) of the distribution of the weighted sum of N assets : 

N 

Ps{s)= / dx P(x)5( V - s) . (60) 

Using the change of variable ([l2|), allowing us to go from the asset returns Xj's to the transformed returns 

Yi's, we get 

Ps{s) = / dy e-'^y'^-'y 5{Y, w,sgn{y,)r\yh - s) . (61) 

Taking its Fourier transform Ps{k) = f dsPs{s)e^'''' , we obtain 

Psik) = ^^^^^ [ dy e-iy^^-V+..E-.^.^.n(..)/-(.?) , (62) 

where Ps is the characteristic function of Ps- 

In the particular case of interest here where the marginal distributions of the variables Xi's are the modified 
Weibull pdf, 

f-Hy^)=X^\^\'' (63) 

with 

qi = 2/ci , (64) 

the equation (62) becomes 

The task in front of us is to evaluate this expression through the determination of the moments and/or 
cumulants. 



5.2.1 Case of independent assets 



In this case, the cumulants can be obtained explicitely []26|]. Indeed, the expression (|65|) can be expressed as 
a product of integrals of the form 



+ 00 



We obtain 



du e 



N 



HkwiXi 



C2n = ^c{n,q.i){Xi 



iWi) , 



\2n 



and 



c(n,gi) = (2n)! 



' (2n-2p)!7ri/2 



p=0 



2!7rV2 



n 



2!7rV2 



(66) 



(67) 



(68) 



Note that the coefficient c(n, qi) is the cumulant of order n of the marginal distribution (|^) with c = 2/qi 
and X = 1- The equation ( |67| ) expresses simply the fact that the cumulants of the sum of independent vari- 
ables is the sum of the cumulants of each variable. The odd-order cumulants are zero due to the symmetry 
of the distributions. 
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5.2.2 Case of dependent assets 



Here, we restrict our exposition to the case of two random variables. The case with N arbitrary can be 
treated in a similar way but involves rather complex formulas. The equation ( |65[ ) reads 

1 



Ps{k) 



dyidy2 exp 



21x^1 V - p 

+X2W2sgn{y2 

and we can show (see appendix ^ that the moments read 
with 

7gig2(2n,2p) 



yi 



V2 



91 



+ 



y2 


92-^- 


V2 





(69) 



(70) 



2p 2(n-p) r {qip + ^) r {q2{n - p) + ^) 
Xi X2 



2F1 -qip, -q2{n-p); p' 

^q^q^lZn, Zp + L) — ZXi X2 1 P2^l[-qiP 



02 + 1 3 2 

,-q2{n-p) + -^—'^2'^ 



, (VI) 

2 ' 
(72) 



where 2F1 is an hypergeometric function. 

These two relations allow us to calculate the moments and cumulants for any possible values of qi = 2/ci 
and q2 = 2/c2- If one of the qi's is an integer, a simplification occurs and the coefficients 'y{n,p) reduce to 
polynomials. In the simpler case where all the qi's are odd integer the expression of moments becomes : 



with 



p=0 



(2n) 
^2p 



n 



P 



min{qip,q2{n~p)} 



(^^1X1)^(^^2X2) 



{2p)\c',;^-^\2{n - p) - 1)U 

0, 
, 



n—p 



s=0 



(2n) 
2p+l 
{2n+l) 
hp 

(2n+l) 



(2n)! 



(73) 



2(2"-2p)(2p)!(n -p)! 



(2p + l)!C72t+7)(2(n-p)-l)!! 



(2n + 1)! 



2(2n-2p)(2j5+ i)!(n-p)! 



(74) 

(75) 
(76) 

(77) 



5.3 Non-symmetric assets 



In the case of asymmetric assets, we have to consider the formula (25 -2^, and using the same notation as in 
the previous section, the moments are again given by ( |70| ) with the coefficient j{n,p) now equal to : 



lin,p) 



(-i)"(xr)'^(x2-)"~^ 



r(0 + iV|^g2-(n-P) + i) 



QiP 12 - P) 1 2V 

' 2 ' 2 '2'^ r 
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+ 



{-inxinxtr-" 



47r 



9i P + 1^ P qiin-p) + 1 



47r 



-2r,^ + 



+ 



47r 



p 2F1 



2F1 



2F1 



2F1 



'2'^ 



2 ' 


2 






(n - p) - 


1_ 




2 


1 




q2in- 


P) 


. 2 ' 


2 




1 ^2" 


(n - p) - 


1_ 




2 




qtp 







qjp-'^ qt{n-p)-l 3 2 
2 ' 2 '2'^ 



+ 



+ 



2^^ ) + 



+ 



'2'^ ) + 



This formula is obtained in the same way as the formulas given in the symmetric case. We retrieve the 



formula ( [71[ ) as it should if the coefficients '+' are equal to the coefficients '- 



5.4 Empirical tests 



Extensive tests have been performed for currencies under the assumption that the distributions of asset 
returns are symmetric [|2^]. 

As an exemple, let us consider the Swiss franc and the Japanese Yen against the US dollar. The calibration 
of the modified WeibuU distribution to the tail of the empirical histogram of daily returns give {qcHF = 
1.75, ccHF = 1-1^,XCHF = 2.13) and {qjpy = 2.50, cjpy = 0.8, xjpy = 1-25) and their correlation 
coefficient is p = 0.43. 



The figure [Ij plots the excess kurtosis of the sum wchfxchf + wjpyxjpy as a function of wchf, de- 
noted w in the figure, with the "no-short" constraint wchf + wjpy = 1- The thick solid line is determined 
empirically, by direct calculation of the kurtosis from the data. The thin solid line is the theoretical predic- 
tion using our theoretical formulas with the empirically determined exponents c and characteristic scales x 
given above. While there is a non-negligible difference, the empirical and theoretical excess kurtosis have 
essentially the same behavior with their minimum reached almost at the same value of wchf- 

Three origins of the discrepancy between theory and empirical data can be invoked. First, as already pointed 
out in the preceding section, the modified Weibull distribution with constant exponent and scale parame- 
ters describes accurately only the tail of the empirical distributions while, for small returns, the empirical 
distributions are close to a Gaussian law. While putting a strong emphasis on large fluctuations, cumulants 
of order 4 are still significantly sensitive to the bulk of the distributions. Moreover, the excess kurtosis is 
normalized by the square second-order cumulant, which is almost exclusively sensitive to the bulk of the 
distribution. Cumulants of higher order should thus be better described by the modified Weibull distribution. 
However, a careful comparison between theory and data would then be hindered by the difficulty in esti- 
mating reliable empirical cumulants of high order. The second possible origin of the discrepancy between 
theory and data is the existence of a weak asymmetry of the empirical distributions, particularly of the Swiss 
franc, which has not been taken into account. The figure also suggests that an error in the determination of 
the exponents c can also contribute to the discrepancy. 



15 



In order to study investigate the sensitivity with respect to the choice of the parameters q and p, we have 
also constructed the dashed line corresponding to the theoretical curve with p = (instead of p = 0.43) and 
the dotted line corresponding to the theoretical curve with qcHF = 2 rather than 1.75. Finally, the dashed- 
dotted line corresponds to the theoretical curve with qcHF = 1.5. We observe that the dashed line remains 
rather close to the thin soUd line while the dotted line departs significantly when w increases. Therefore, the 
most sensitive parameter is q, which is natural because it controls directly the extend of the fat tail of the 
distributions. 

In order to account for the effect of asymmetry, we have plotted the fourth cumulant of a portfolio compound 



of Swiss Franc and British Pound. On figure [15|, the solid line represents the empirical cumulant while the 
dashed line shows the theorical cumulant. The agreement between the two curves is better than under the 
symmetric asumption. Note once again that an accurate determination of the parameters is the key point 
to obtain a good agreement between empirical data and theoretical prediction. As we can see in the figure 



15, the paramaters of the Swiss Franc seem well adjusted since the theoretical and empirical cumulants are 
both very close when w ~ 1, i.e., when the Swiss Franc is almost the sole asset in the portfolio, while when 
~ 0, the theorical cumulant is far from the empirical one, i.e., the parameters of the Bristish Pound are 
not sufficiently well-adjusted. 



6 Portfolio optimization 

Up to now, we have calculated several measures of risk: the asymptotic distribution of losses and the corre- 
sponding Value-at-Risk to quantify the large fluctuations and the cumulants to evaluate smaller fluctuations. 
We now present different strategies of portfolio allocation constructed using these risk measures. 



6.1 Optimization with respect to the Value-at-Risk 

As already said, in the case of assets with different exponent c, there is no possible diversification of large 
risks. Therefore, the strategy minimizing large risks consists in holding as little as possible of the more 
fat-tailed assets, compatible with the constraint on the average return of the portfolio. Obviously, if only 
large risk control matters, one should only hold the asset with the thinnest tail. 

The more interesting case occurs when the exponents c are the same for the different assets of the portfolio. 
According to the tables |] and ^ this is indeed an relevant situation, when taking into account the error bars 
in the exponents. It is then easy to show, from the results given in section ^ that the VaR is an increasing 
function of the scale parameter x where: 

^ " ..-1 c/2 c/2 '^'^^^ 

when c > 1 and 

= ma-K {wlxi,- ■ ■ ,wlix%] (80) 

i&{l,..N} 

if c < 1. 

Thus, the optimization program to fulfill is : 

X* = infj^j VaR({Ti;i}) 

EiW^i = l (81) 

Ei wi^J'i = At 

where /ij denotes the average return of asset i and p the average return of the portfolio. 
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Case c > 1 : In the case, we can give the analytical form of the efficient frontier, i.e., /i versus x*'^ (^t 
least) in the independent case. Let Ai and A2 denote two Langrange multipliers. We have to solve 



d 



X 



Xi(^Wi-l^ -X2(^Wifi}j^ =0, (82) 



dwk 

which leads to 

( I^KXi)^ ] xfwl^^-' = Ai + X2i2k ■ (83) 



c-1 



Multiplying by zu^ and summing over k, we obtain : 

Ai + Aa/i = X . (84) 

Then, expressing wl from eq.(|^) and accounting for the constraint, we obtain an analytical parametric 
equation of the efficient frontier : 

( (Ai+A2At»)''~^ _ ^* 

E, = ^x* (85) 

\ Ai + A2/u = x* 

Varying Ai and A2, the efficient frontier is delineated. In the Gaussian case c = 2, we retrieve the standard 
Markovitz efficient frontier. 

In the case of correlated assets, the shape of the frontier remains essentially the same but its determination 
requires numerical calculation. 

Case c < 1 : In this case, the general minimization problem must be again solve numerically. As an 
example, we consider the simple case where the portfolio is made of only two assets. This case is analyticaly 
soluble. Indeed, the program 

inft„*G[o,i]{suPig{i,..7V}{wf X?, • • • , wlixl;}} 

Ei"^^ = l (86) 

Ei wi^J'i = iJ' 

simplifies into : 

2 "1 





-/^2 \ 2 








Xi, 








-/^2y 







inf <'sup<;(^^^) Xl\\- (87) 

The thick line on figure |l^ represents the efficient frontier while the dashed lines represent fi = f{xi) ^^'^ 

1^ = /(X2)- 



6.2 Optimization with respect to the cumulants 

A priori, all order cumulants can be considered and each of them embodies a certain measure of risk. 
Obviously, the larger the cumulant order, the more sensitive is the cumulant to large fluctuations. We can 
expect risk control based on the high-order cumulants will be equivalent to risk-control based on the VaR 



approach [26|. In this section, we will only focus on relatively low-order cumulants, because they are 
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expected to add valuable information on top of the risk structure already provided by the VaR approach. 
By lower-order, we mean cumulants of order n = 2 to n = 6. To justify this approach, we have evaluated 
the cumulants for several portfolios and it is very clear that, in most cases for n > 6, they become strictly 
monotonous functions of the asset weights within the interval [0, 1]^. Therefore, the minimization of risks 
which would correspond to minimizing cumulants of order larger than n = 6 would lead once more to hold 
as Uttle as possible of the more fat-tailed assets. 

The portfolio optimization along the lines of [^, 25, 2|] corresponds to minimizing a given cumulant C„ of 
order n of the distribution of returns of the portfoUo in the presence of a constraint on the return as well as 
the "no-short" constraint: 

infu,,e[o,i] Cn{{wi}) 

EiWi = '^ (88) 

In the most general case, we have to use extentions for > 2 of the formulae given in section ^ and perform 
a numerical optimization. Interesting results can be observed in simpler situations, as we now obtain. 

Figure |l^ and 18 show the generalized efficient frontiers using C2 (Markovitz case), C4 or Cq as relevant 
measures of risks, for two portfolios composed of two stocks : IBM and Hewlett-Packard in the first case 
and IBM and Coca-Cola in the second case. 

Obviously, given a certain amount of risk, the mean return of the portfolio changes when the cumulant 



considered changes. It is interesting to note that, in figure 17, the minimisation of larger risks, i.e., with 



respect to Cq, increases the average return while, in figure 18 , the minimisation of larger risk lead to decrease 
the average return. 

This allows us to make precise and quantitative the previously reported empirical observation that it is 
possible to "have your cake and eat it too" We can indeed give a general criterion to determine under 
which values of the parameters (exponents c and characteristic scales x of the distributions of the asset 
returns) the average return of the portfolio may increase while the large risks decrease at the same time, 
thus allowing one to gain on both account (of course, the small risks quantified by the variance will then 
increase). For two independent assets, assuming that the cumulants of order n and n + k oi the portfoUo 
admit a minimum in the interval ]0, 1[, we can show that 

< < <+fc (89) 

if and only if 



fii - fi2 > and 



Cn(l) 
Cn{2) 



1 

n-l 



> 



Cn+fc(2) 



(90) 



fJ-i — < and 



Cn{l)\- 



'1 ^ ,^C7„+fc(l)V+^- 



Cn+k{2) 



(91) 



where /x* denotes the return of the portfolio evaluated with respect to the minimum of the cumulant of order 
n and C„(i) is the cumulant of order n for the asset i. 

The proof of this result and its generalisation to > 2 are given in appendix ^ In fact, we have observed 
that when the exponent c of the assets remains sufficiently different, this result still holds in presence of 
dependence between assets. 
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For the assets considered above, we have found fxiBM = 0.13, hhwp = 0.07, fXKO = 0.05 and 



C2{IBM) 
C2{HWP) 

C2{IBM) 
C2{KO) 



1 7R ^ / C4{IBM) \ 3 



0.96 < 



C4iIBM) \ 3 
CiiKO) ) 



1.03 > 



1.01 < 



(C^{IBM)_y 

C&{IBM)\^ 
CeiKO) J 



0.89 



1.06 



(92) 
(93) 



which shows that, for the portfolio IBM / Hewlett-Packard, the efficient return is an increasing function of 
the order of the cumulants while, for the portfolio IBM / Coca-Cola, the inverse phenomeon occurs. This is 



exactly what is shown on figures 17 and Mq. 



The underlying intuitive mechanism is the following: if a portfolio contains an asset with a rather fat tail 
(many "large" risks) but narrow waist (few "small" risks) with very little return to gain from it, minimizing 
the variance C2 of the return portfolio will overweight this asset which is wrongly perceived as having little 
risk due to its small variance (small waist). In contrast, controlling for the larger risks quantified by C4 or 
Cq leads to decrease the weighing of this asset on the portfoio, and correspondingly to increase the weight 
of the more profitable assets. We thus see that the effect of "both decreasing large risks and increasing 
profit" appears when the asset(s) with the fatter tails, and therefore the narrower central part, has(ve) the 
smaller overall return(s). A mean-variance approach will weight them more than deemed appropriate from 
a prudential consideration of large risks and consideration of profits. 



7 Conclusion 

We have presented new analytical results on and empirical tests of a general framework for a portfolio 
theory of non-Gaussian risks with non-linear correlations. We have shown that the concept of efficient 
frontiers can be generalized by using other measures of risks than the variance (cumulant of order 2), for 
instance the value-at-risk and the cumulants of order larger than 2. This work opens several novel interesting 
avenues for research. One consists in extending the Gaussian copula assumption, for instance by using the 
maximum-entropy principle with non-extensive Tsallis entropies, known to be the correct mathematical 
information-theoretical representation of power laws. A second line of research would be to extend the 



present framework to encompass simultaneously different time scales r in the spirit of [ |19| ] in the case of a 
cascade model of volatilities. 
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A Asymptotic expansion of the wealth distribution for a portfolio made of 
assets with same exponent c 



A.l Case of independent assets 

We start with the definition S = YliLi '^i^i the corresponding equation for its probability density 
function : 

Ps{S)o^ J dxi---dxN e'^^^^^^'^^'^U (^S -Y^Wix}j , (94) 
where we have used the parameterization 

Piixi) = Ci e-^'("') , (95) 

and 



h{xi)= (^) , withol. (96) 



All equations in (Q) are from — oo to +oo. The delta function expresses the constraint on the sum. We need 
the following conditions on the functions fi 

(i) fi{xi) — > +00 sufficiently fast to ensure the normalization of the pdf's. 

(ii) fl'{xi) > (convexity), where /" denotes the second derivative of /. 

(iii) liuix^oox'^fix) = +00. 

Under these assumptions, the leading order expansion of PsiS) for large S and finite > 1 is obtained by 
a generalization of Laplace's method which here amounts to remark that the set of x*'s that maximize the 
integrand in (Q) are solution of 

flix*) = a{S)wi , (97) 

where a{S) is nothing but a Lagrange multiplier introduced to minimize the expression J2iLi f-ii^i) under 
the constraint X^^j^ WiXi = S. 

Expanding fi{xi) around x* yields 

nixi) = fi{xi) + fi{xi)hi + ^/f + \h{xr)^'^hi + ^jl'\x*)ht ■ ■ ■ m 

where the hi = Xi — x* obey the condition 

N 



Y,w^h, = 0. (99) 



i=l 

Taking into account equation (p7|), we obtain 

_ (100) 

and with 1 



E M^^) = E /^«) + ^ E + ^ E /f «)^' + ^ E + • • • (101) 
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Thus exp(— J2 fi{xi)) can be rewritten, up to the /i^-order, as follows : 



(102) 



Then expression ( 94 ) becomes 



Using the fact that 



2tt 



(103) 



(104) 



we obtain 



P5(5)oce-S" 1^ jdh,---dhN e-h^Uif^i^m+^^f^-'.'^^) 



(105) 



and 

^5(5) oc e- 



j dhi ■ ■ ■ dh^ e 



(106) 



Let us denote by the average over hi with respect to a Gaussian distribution whose mean is —ikj^^^ 



and variance : 1 / yjffixf) and by the average over with respect to a Gaussian distribution whose mean 



is and variance : 1/a f"(x*) ■ ^^^^ these notations, we can rewrite the equation above as follows 



Ps{S) oc 



(2vr)- 



24 



(107) 



The term (/if)? (and every odd order term in hi) only involves odd order term in k and thus {{h^)i)k vanishes, 
so a simplification occurs : 



Ps{S) oc 



(2vr) 



JV-l 
2 



^^=1 WW) 



i + iE/i'^«)((^')^)^ + --- 



24 



(108) 



The evaluation of {{hj)i)i. leads to 



+ 



3wf 



(109) 
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Using equation (37), we obtain 



-WiXi 



1 

c-l 



which, together with the constraint J2iLi '^i^t = "S*. gives 

I 



where 



N 



-s, 



c ^ — ^ c 

1=1 



(110) 



(111) 



(112) 



Thus, X* oc S and ^(x*) « 5'=-" So {{hf)i)i, a 5^-2'= and ^ Hx*)((/i4),)fc ^ 5-c_ 
Let us now analyze the impact of the higther order terms. First of all, consider the G*'' order term 

where we have omitted the numerical constants. More generaly, for the m*'^ order term, we will obtain a 
contribution of the form 

i 

and many other contributions of the form 
with Pi + ■ ■ ■ + pn = m and pi > 3. 

Thus, as a first step, we have to evaluate the n*'* moment (/i") of a Gaussian variable h. Let us denote m 
and a the mean and variance of this variable. Using the well known result 



(113) 



exp (-t^ + 2tx) = V — , 
where Hn{-) is the Hermite polynomial of order n, we get 



im 



and so 



3 jn 



Wjk 



(116) 



(117) 



(118) 



As a second step, let us evaluate the term J2i fi'^"^\xi){{^i"^)'i)k (recall that odd order terms vanish). We 
have 



{{hr)j)k = i-ir I dk 
(-ly 



E 



TT — / \ m 1 I 



Wjk 



2m 



\ 



(119) 
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The integral can be calculated exactly but this is not useful as the only interesting thing is that the argument 
in the polynomial is independent of S. Indeed, the terms //'(x,*)'s are proportional to S"*^"^ and cancel 
out in the argument of H2m- Thus the integral behaves like a constant, and the dependence with respect to 
S in {{h'j"^)j)k is only given by the prefactor f"{x*)~'^ : 

oc -^^^ . (120) 

As ff'^\x^) behaves like S"^"^"*, we obtain 

i 

where m > 2. Therefore, these terms decrease at least as fast as (with c > 1). 
The third and final step consists in evaluating the "mixed"-terms. We have 

((^Dii • • • (^Di )k oc 7^ ^(-ffj,! ■■■Hr, )k, (122) 

" " (2//:Kr/^---(2//:«)r'' 

where, as in the second step, the average of the Hp^'s over k is independent of S. Thus, 

■ ■ ■ {h^hJk oc oc . (123) 

Taking into account the fact that 

fii'^ ■ ■ ■ AT^ °^ ^"'^^ ■ ■ ■ ^"'^^ °^ 5"^-2"* , (124) 

we are led to 

E ft\<) ■ ■ ■ ft\<mih. ■ ■ ■ mi.h oc . (125) 

As H \-Pn = 2m and pi > 3, we have 3n < 2m. Som — n>^>l and the "mixed"-terms decrease 

at least like S~'^. 

Finally, we can conclude that terms of order higher than 1 are sub-dominant for large S which provides our 
exact asymptotic result: 

Ps{S)oc . ^^"^"^ e-^l.M^:)(i + o(l\) . (126) 



^From this, we obtain 



= ^^^5^ ^ E/.(x:) = (^-]\ (127) 



X 



f -,^ 2 ^ ( ( ^\ \^ N 

= =^ n /no = %4s'^'' n -t'^' (128) 
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Finally, we find : 



iN 



xi 



N 



X 



2-c 



27r 



c(c-l) 



2-c 



N-l 
2 



(130) 



Note that this calculation does not require that the functions /j's should be symmetric. Therefore, the result 
still holds for asymmetric functions. 

More over, it is easy to show that when we consider pdf's whose equations are given by gi{x) e~^'^^^ the 
same result holds as long as lngi{x) remains sub-dominant with respect to fi{x). Therefore, we just have 
to multiply equation (|130D by Y\ Qiix*) to obtain the corresponding correct result. Thus, for the pdf's given 



by (|24|), the equation ( |130| ) becomes : 



P{S) 



iV+l 
C 2 



N 



iV+l 

X2 2 (c 



jWiwiXi 

1=1 



c-l 
\ 2 



2{c-l) 



\s\ 

X 



c-2 



•|S|^ 



(131) 



A.2 Case of dependent assets 



We assume that the marginal distributions are given by the modified WeibuU distributions: 



Pi{Xi) 



_i c_ 



(132) 



Under the Gaussian copula assumption, we obtain the following form for the mutivariate ditribution 



• • ■ ,Xn) 



N c/2-1 
X- 



2JV^iV/2^1_i ^c/2 



n 



■ exp 



c/2 / ^ c/2 



Xi 



Xj 



(133) 



Let 



(134) 



We have to maximize / under the constraint ^ WiXi = s. As for the independent case, we introduce a 
Lagrange multiplier A which leads to 



(135) 



and using the constraint "^WiX* = S : 



(136) 
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Inspired by the results obtained for independent assets and by dimensional arguments, we assert that x* is 
proportional to S and A proportional to S'^~^ : 



A = KS"-^ , 

where ai and K depend on {wj} and {xj} but are independent of S. 
Therefore, 



c/2 c/2 ^ 



and 



Thus 



(137) 
(138) 



(139) 
(140) 



(9/ 1 / 



E^.-l c/2 c/2^c , 1 , , , 



where, as in the previous section, hi = Xi — x*. 



(141) 
(142) 



It is easy to check that the n*^-order derivative of / with respect to the Xj's evaluated at {x^} is proportional 
to /S*^"". In the sequel, we will use the following notation : 



dxi^ • • • <9xj„ 



K} 



n—in 



(143) 



We can write : 



/(xi, . . . , x^) = lA^Vf af 5= + ^ E ^if^^h^ + ^ E ^lfkh^h,h, + 

j,k ij ijk 



(144) 



up to the fourth order. This leads to 

Ei,-1 c/2 c/2 „„ f 



dhi---dhNe — ^^i^A^i 



(2) 



S' 



c-3 



1 + VE^S^^^^-^^ + 



6 



(145) 



Using the same method as in the previous section : 



Mlfhihj-ikY^jWjkj ^ 
(3), 



i + ^E^S^^^A + --- 



ijk 



(146) 
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or in vectorial notation : 



P{S) oc e 



ETT—l c/2 c/'2 



dk 
2^ 



dh e 2 



(147) 



Let us peri'orm the following standard change of variables : 



h = h'- 



Sc-2 



w , 



(M^^) ""^ exists since / is assumed convex and thus M^^^ positive) : 



c-2 



-h*M(2)h + zA;w*h 



h'*M(2)h' + -^w*M(2)-V 
2 25^-2 



This yields 



/ 



X / dhe 2 



^c-3 



e 2SC- 



J^w*M(2)-V 



(148) 



(149) 



ijk 



. (150) 



Denoting by {■)h the average with respect to the Gaussian disttibution of h and by {■)k the average with 
respect to the Gaussian distribution of fc, we have : 



P{S) oc . h^:^^^ (2.5-^)^ ^'^^..y^'^'^'s^ ^ 

^ V wtM(2)-V 



ijk 



ijkl 



(151) 



We now invoke Wick's theorem, which states that each term {{hi • • • hp)Yi)}^ can be expressed as a product 
of pairwise correlation coefficients. Evaluating the average with respect to the symmetric distribution of k, 
it is obvious that odd-order terms will vanish and that the count of powers of S involved in each even-order 
term done for independent assets remains true. So, up to the leading order : 



P{S) oc \l _ (27r5^~^)~ e "^j^kVjk ''j ^ 



w 



tM(2)" 



(152) 



w 
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B Calculation of the moments of the distribution of portfolio returns 



Let us start with equation (pS) in the 2-asset case 



Psik) 



1 



+ ik (^Xiwisgn{yi) 



dyidy2 exp 

+X2W2sgn{y2 

Expanding the exponential and using the definition of moments, we get 



yi 



V2 



91 



+ 



y2 


92-^ - 


V2 





(153) 



27ryr^ 



== / dyidy2y2\ 1X1X2 ^w^w^ '^sgn{yiY 

xsgn{y2T~P 



yi 



V2 



HP 



V2 



q2{n-p) 



(154) 



Posing 

79192 

this leads to 



27rVl - 



dyidy2 sgn{yiY 



yi 



V2 



91P 



sgn{y2T-P 



y2 



V2 



92(n-p) 



^-lyty-iy 



= ( J ^1^2 ^79192 

p=0 ^■'^^ 



(n,p) . 



(155) 



(156) 



Let us defined the auxiliary variables a and /3 such that 



a 



)12 '^"^ 



-iV-^)l2 = -iV-')21 - 



(157) 



Performing a simple change of variable in ( |155| ), we can transform the integration such that it is defined 
solely within the first quadrant (yi > 0, y2 > 0), namely 



79192 (^'P) = Xl^X2' 



n r+00 



^ , ..... J^y7^y'""^-t(^?-^^) X 



-p i + (-ir 



This equation imposes that the coefficients 7 vanish if n is odd. This leads to the vanishing of the moments 
of odd orders, as expected for a symmetric distribution. Then, we expand e^^^^^ _)_ (_i)Pg-/3yij/2 series. 
Permuting the sum sign and the integral allows us to decouple the integrations over the two variables yi and 

y2- 



1 -I- f— 1 

79192 = XlV~'--j=l= + (-IF^I^ 



+00 „,9iP+s ^ \ 

2— / 



<j2(n-p)+s ^ 

dyo ^ e-t^? 

2 2 



(159) 
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This brings us back to the problem of calculating the same type of integrals as in the uncorrelated case. 
Using the expressions of a and /?, and taking into account the parity of n and -p, we obtain: 



n _ rt2\gip+g2(n-p)+| +2? C9o^2s / -i \ 



xr( g2(n-p) + s + M , (160) 



7^,52 (2n,2p + l) - XI X2 ^ (2g + l)! 



1 

X 

('2.S -I- I H 

s=0 



xr(gip + s + l + |-)r(g2(n-p) + s + l-|) . (161) 

Using the definition of the hypergeometric functions 2^1 l[l]l, and the relation (9.131) of we finally 
obtain 

TT 

X2Fi(-5ip-^,-52(n-p) + ^;^;p^) . (163) 
In the asymmetric case, a similar calculation follows, with the sole difference that the results involves four 



terms in the integral (158) instead of two. 
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C Conditions under which it is possible to increase the return and decrease 
large risks simultaneously 



We consider N independent assets {1 • • • N}, whose returns are denoted by /xi • • • fijy. We aggregate these 
assets in a portfolio. Let wi - ■ ■ wn be their weights. We consider that short positions are forbiden and that 
J2- Wi = 1. The return fi of the portfolio is 

N 

lx = ^Wiiii. (164) 

i=l 

The risk of the portfolio is quantized by the cumulants of the distribution of /x. 

Let us denote /x* the return of the portfoUo evaluated for asset weights which minimize the cumulant of 
order n. 

C.l Case of two assets 

Let Cn be the cumulant of order n for the portfoho. The assets being independent, we have 

Cn = C„(l)w;i" + C„(2)ix;2", (165) 
= C„(l)w;i" + C„(2)(l-ii;i)". (166) 

In the following, we will drop the subscript 1 in ix;!, and only write w. Let us evaluate the value ix; = xx;* at 
the minimum of Cn, n > 2 : 



Cn(2) V 

and assuming that C„(l)/C„(2) > 0, we obtain 



^ C„(l)xx;"-i-C„(2)(l-xx;)"-i =0, (167) 
Cnil)_(l-w*r-^ (168) 



xx;* = ^-(^^^ . ■ (169) 

Cn{l)—^+Cn{2)—^ 

This leads to the following expression for /x* : 

, _ /xi-C„(2)^+/X2-Cn(l)^ 

Mn — 11- 

C„(l)^+C„(2)^ 
Thus, after simple algebraic manipulations, we find 

Hl < ^ - 112) (c„(l)^C„+fc(2)^ - C„(2)^C„+fc(l)^) > 0, (171) 
which concludes the proof of the result annouced in the main body of the text. 
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C.2 General case 



We consider a portfolio with N independent assets. Assuming that the cumulants C„(i) have the same sign 
for all i, the minimum of C„ is obtained for a portfolio whose weights are given by 

nf^.gn(j)^ 

= ^> (172) 



and we have 



< = _ , . 1 ■ (173) 



Indeed, the following conditions hold: 



N 



Y,Cnii)w'l-'dxi = 0, (174) 

N 

dxi = 0. (175) 



i=l 

N 



i=l 



Introducing a Lagrange multiplier A, we obtain 



Assuming that all the C„ (i) have the same sign, we can find a A such that all the Wi are real and positive : 



Wi = —. (177) 

If some Cn{i) are positive and others negative, the set of assets we consider is not compatible with a global 
mmimum on ]0, Then, we have to spht the portfolio into two sub-portfolios constituted of assets whose 
cumulants have the same sign and perform the minimization of the corresponding cumulant of each subset. 

There is not simple condition that ensure /U* < /i *^^. The simplest way to compare /^* and ^*^^ is to 



calculate diretly these quantities using the formula ( |173| ). 



30 



References 



Abramovitz, E. and I.A. Stegun, Handbook of Mathematical functions (Dover Publications, New 
York 1972). 

Andersen, J.V., and D. Sornette, Have your cake and eat it too: increasing returns while lowering 



large risks! preprint at http://xxx.lanl.gov/abs/cond-mat/9907217 



Artzner, P., F. Delbaen, J.M. Eber and D. Heath, Coherent measures of risk. Math. Finance 9, 203-228 
(1999). 

Bouchaud, J.-P, D. Sornette, C. Walter and J.-P. Aguilar, Taming large events: Optimal portfolio 
theory for strongly fluctuating assets. International Journal of Theoretical and Applied Finance 1, 
25-41 (1998). 

Embrechts, P., C. Kluppelberg and T. Mikosh, Modelling extremal events (Springel-Verlag, Applica- 
tions of Mathematics 33, 1997). 

Embrecht, P., A. McNeil and D. Straumann, Correlation and Dependence in risk management: prop- 
erties and pitfalls. Proceedings of the Risk Management workshop, October 3, 1998 at the Newton 
Institute Cambridge, Cambridge University Press. 

Frisch, U. and D. Sornette, Extreme Deviations and Applications, J. Phys. I France 7, 1155-1171 
(1997). 

Gopikrishnan, P., M. Meyer, L.A. Nunes Amaral and H.E. Stanley, Inverse cubic law for the distri- 
bution of stock price variations, European Physical Journal B 3, 139-140 (1998). 

Gradshteyn, I.S. and I. M. Ryzhik, Table of integrals Series and Products, Academic Press, 1965. 

Hill, B.M., A Simple General Approach to Inference about the Tail of a Distribution, Annals of 
statistics, 3(5), 1163-1174 (1975). 

Jorion, P., Value-at-Risk: The New Benchmark for Controlling Derivatives Risk (Irwin Publishing, 
Chicago, IL, 1997). 

Karlen, D., Using projection and correlation to approximate probability distributions. Computer in 
Physics 12, 380-384 (1998). 

Laherrere, J. and D. Sornette, Stretched exponential distributions in nature and economy : "fat tails" 
with characteristic scales, European Physical Journal B 2, 525-539 (1998). 

Lindskog, F., Modelling Dependence with Copulas, 

http : // www.risklab.ch/ Papers.html^MTLindskog 



Litterman, R. and K. Winkehnann, Estimating covariance matrices (Risk Management Series, Gold- 
man Sachs, 1998). 

Lux, T., The stable Paretian hypothesis and the frequency of large returns: an examination of major 
German stocks. Applied Financial Economics 6, 463-475 (1996). 

Markovitz, H., Portfolio selection : Efficient diversification of investments (John Wiley and Sons, 
New York, 1959). 

Merton, R.C., Continuous-time finance, (Blackwell, Cambridge, 1990). 



31 



[19] Muzy, J.-R, D. Sornette, J. Delour and A. Arneodo, Multifractal returns and Hierarchical Portfolio 
Theory, Quantitative Finance 1 (1), 131-148 (2001). 

[20] Nelsen, R.B. An Introduction to Copulas. Lectures Notes in statistic, 139, Springer Verlag, New York 
(1998). 

[21] Pickhands, J., Statistical Inference Using Extreme Order Statitstics, Annals of Statistics, 3(1), 119- 
131 (1975). 

[22] Rao, C.R., Linear statistical inference and its applications, 2d ed. (New York Willey, 1973). 

[23] Sornette, D., Large deviations and portfolio optimization, Physica A 256, 251-283, 1998. 

[24] Sornette, D., Critical Phenomena in Natural Sciences, Chaos, Fractals, Self-organization and Dis- 
order: Concepts and Tools, 432 pp., 87 figs., 4 tabs (Springer Series in Synergetics, August 2000) 
ISBN 3-540-67462-4 

[25] Sornette, D., J. V. Andersen and P. Simonetti, Portfolio Theory for "Fat Tails", International Journal 
of Theoretical and Applied Finance 3 (3), 523-535 (2000) 

[26] Sornette, D., P. Simonetti, J.V. Andersen, (/)^-field theory for portfolio optimization : "fat-taUs" and 
non-linear correlations. Physics Reports 335 (2), 19-92 (2000). 

[27] TsalUs, C. Possible generahzation of Boltzmann-Gibbs statistics, J. Stat. Phys. 52, 479-487 (1988); 



for updated bibliography on this subject, see fittp : //tsallis.cat.cbp f.br/biblio.htm 



32 





< X+ > 


< c+ > 


X+ 




< X- > 


< c_ > 


X- 


c_ 


CHF 


2.45 


1.61 


2.33 


1.26 


2.34 


1.53 


1.72 


0.93 


DEM 


2.09 


1.65 


1.74 


1.03 


2.01 


1.58 


1.45 


0.91 


JPY 


2.10 


1.28 


1.30 


0.76 


1.89 


1.47 


0.99 


0.76 


MAL 


1.00 


1.22 


1.25 


0.41 


1.01 


1.25 


0.44 


0.48 


POL 


1.55 


1.02 


1.30 


0.73 


1.60 


2.13 


1.25 


0.62 


THA 


0.78 


0.75 


0.75 


0.54 


0.82 


0.73 


0.30 


0.38 


UK? 


1.89 


1.52 


1.38 


0.92 


2.00 


1.41 


1.82 


1.09 



Table 1: Table of the exponents c and the scale parameters x for different currencies. The subscript "+" 
or "-" denotes the positive or negative part of the distribution of returns and the terms between brackets 
refer to parameters estimated in the bulk of the distribution while naked parameters refer to the tails of the 
distribution. 
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^ C-j- ^ 
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< X- > 




X- 




A A/TAT 


IZ.^ / 


1 .oZ 


8 7'^ 
o. / J 


n QQ 


1 1 OA 


l.uu 


8 11 


n Q8 


p 

V 


0. 


1 70 


9 'XA 
Z. 


u.uu 


^ QQ 

J.yy 


1 70 
1. /U 


n AC\ 


U. JU 


FA/TP 




1 

1 .Oj 




1 .J J 


1 1 AA 


1 




n '^7 
U.J / 


r 


7 '^7 


1 . JZ 




1 .OH" 


f\ flQ 

\j,\Jy 


1 Q1 


Q7 

J .y 1 


1 ^4 


GE 


5.21 


1.89 


1.81 


1.28 


4.80 


1.81 


4.31 


1.16 


GM 


5.78 


1.71 


0.63 


0.48 


5.32 


1.89 


2.80 


0.79 


HWP 


7.51 


1.93 


4.20 


0.84 


1.26 


1.76 


1.66 


0.52 


IBM 


5.46 


1.71 


3.85 


0.87 


5.07 


1.90 


0.18 


0.33 


INTC 


8.93 


2.31 


2.79 


0.64 


9.14 


1.60 


3.56 


0.62 


KO 


5.38 


1.88 


4.46 


1.04 


5.06 


1.74 


2.98 


0.78 


LU 


10.46 


2.02 


7.12 


0.98 


10.16 


1.79 


2.11 


0.43 


MDT 


6.82 


1.95 


6.09 


1.11 


6.49 


1.54 


2.55 


0.67 


MRK 


5.36 


1.91 


4.56 


1.16 


5.00 


1.73 


1.32 


0.59 


MSFT 


8.14 


2.19 


2.11 


0.58 


1.11 


1.60 


0.67 


0.38 


PEE 


6.41 


2.01 


5.84 


1.27 


6.04 


1.70 


0.26 


0.35 


PG 


4.86 


1.83 


3.53 


0.96 


4.55 


1.74 


2.96 


0.82 


QCOM 


15.15 


1.70 


14.76 


1.40 


13.02 


1.86 


10.17 


1.07 


SBC 


5.21 


1.97 


1.26 


0.59 


4.89 


1.59 


1.56 


0.60 


SUNW 


11.54 


1.94 


6.91 


0.90 


10.96 


1.70 


5.89 


0.77 


T 


5.13 


1.48 


2.86 


0.76 


4.60 


1.70 


1.87 


0.56 


TXN 


9.06 


1.78 


4.07 


0.72 


8.24 


1.84 


2.18 


0.54 


WCOM 


9.80 


1.74 


11.01 


1.56 


9.09 


1.56 


2.86 


0.58 


WMT 


7.41 


1.83 


5.81 


1.01 


6.80 


1.64 


3.75 


0.78 



Table 2: Table of the exponents c and the scale parameters x for different stocks. The subscript "+" or "-" 
denotes the positive or negative part of the distribution and the terms between brackets refer to parameters 
estimated in the bulk of the distribution while naked parameters refer to the tails of the distribution. 
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Bimodal Example PDFs 




f(x) I I 

Figure 1: Schematic representation of the nonlinear mapping Y = u{X) that allows one to transform a 
variable X with an arbitrary distribution into a variable Y with a Gaussian distribution. The probability den- 
sities for X and Y are plotted outside their respective axes. Consistent with the conservation of probability, 
the shaded regions have equal area. This conservation of probability determines the nonlinear mapping. 
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Figure 2: Quantile of the normalized sum of the Gaussianized returns of the Swiss Franc and The British 
Pound versus the quantile of the Normal distribution, for the time interval from Jan. 1971 to Oct. 1998. 
Different weights in the sum give similar results. 
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KO-PG 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Cumulative Normale Ditribution 



Figure 3: Quantile of the normalized sum of the Gaussianized returns of Coca-Cola and Procter&Gamble 
versus the quantile of the Normal distribution, for the time interval from Jan. 1970 to Dec. 2000. Different 
weights in the sum give similar results. 
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MRK-GE 




0.4 0.5 0.6 
Cumulative Normal Ditribution 



Figure 4: Quantile of the normalized sum of the Gaussianized returns of Merk and General Electric versus 
the quantile of the Normal distribution, for the time interval from Jan. 1970 to Dec. 2000. Different weights 
in the sum give similar results. 
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Figure 5: Cumulative distribution of = y*V~^y versus the cumulative distribution of chi-square (de- 
noted x^) with two degrees of freedom for the couple Swiss Franc / British Pound, for the time interval from 
Jan. 1971 to Oct. 1998. This should not be confused with the characteristic scale used in the definition 
of the modified WeibuU distributions. 
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Figure 6: Cumulative distribution of = y*V~^y versus versus the cumulative distribution of the chi- 
square with two degrees of freedom for the couple Coca-Cola / Procter&Gamble, for the time interval 
from Jan. 1970 to Dec. 2000. This should not be confused with the characteristic scale used in the 
definition of the modified WeibuU distributions. 
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MRK-GE 




Figure 7: Cumulative distribution of = y*V~^y versus versus the cumulative distribution of the chi- 
square with two degrees of freedom for the couple Merk / General Electric, for the time interval from 
Jan. 1970 to Dec. 2000. This should not be confused with the characteristic scale used in the definition 
of the modified WeibuU distributions. 
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MAL + 




Figure 8: Graph of Gaussianized Malaysian Ringgit returns versus Malaysian Ringgit returns, for the time 
interval from Jan. 1971 to Oct. 1998. The upper graph gives the positive tail and the lower one the negative 

tail. The two straight lines represent the curves y = y/2 (tt^) and y = \/2 (-^) 
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UKP + 




UKP- 




Figure 9: Graph of Gaussianized British Pound returns versus British Pound returns, for the time interval 
from Jan. 1971 to Oct. 1998. The upper graph gives the positive tail and the lower one the negative tail. The 

two straight Unes represent the curves y = \/2 f tt^^) and y = \/2 (r^) 
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GE + 




Figure 10: Graph of Gaussianized General Electric returns versus General Electric returns, for the time 
interval from Jan. 1970 to Dec. 2000. The upper graph gives the positive tail and the lower one the negative 

tail. The two straight lines represent the curves y = y/2 (tt^) and y = \/2 (-^) 
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IBM + 




Figure 11: Graph of Gaussianized IBM returns versus IBM returns, for the time interval from Jan. 1970 to 
Dec. 2000. The upper graph gives the positive tail and the lower one the negative tail. The two straight lines 

represent the curves y = \/2 (j^') and y = \/2 (-^) 
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Figure 12: Graph of Gaussianized Wall Mart returns versus Wall Mart returns, for the time interval from 
Sep. 1972 to Dec. 2000. The upper graph gives the positive tail and the lower one the negative tail. The two 

straight lines represent the curves y = -\/2 f 7-^) and y = \/2 (-^) 
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a 



Figure 13: Logarithm In |^7^ j of the ratio of the recurrence time T to a reference time Tq for the recurrence 

of a given loss VaR as a function of a defined by a = VaR* (resp. VaR) is the Value-at-Risk over 

a time interval Tq (resp. T). 
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Excess Kurtosis for CHF / JPY 




CHF 



Figure 14: Excess kurtosis of the distribution of the price variation wchfxchf + wjpy^jpy of the 
portfolio made of a fraction wchf of Swiss franc and a fraction wjpy = 1 — wchf of the Japanese 
Yen against the US dollar, as a function of wchf, denoted w in the figure. Thick soUd line : empirical 
curve, thin soUd Une : theoretical curve, dashed line : theoretical curve with p = (instead of p = 0.43), 
dotted line: theoretical curve with qcHF = 2 rather than 1.75 and dashed-dotted line: theoretical curve with 
qcHF = 1-5. The excess kurtosis has been evaluated for the time interval from Jan. 1971 to Oct. 1998. 
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Fourth Cumulant for CHF / UKP 
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Figure 15: Fourth cumulant for a portfolio made of a fraction wchf of Swiss Franc and 1 — wchf of 
British Pound. The thick soUd Une represents the empirical cumulant while the dotted Une represents the 
theoretical cumulant under the symmetric assumption. The dashed line shows the theoretical cumulant when 
the shght asymmetry of the assets has been taken into account. This cumulant has been evaluated for the 
time interval from Jan. 1971 to Oct. 1998. 
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1 



Figure 16: Efficient frontier (thick line) in the case of two assets with the same exponent lower than 1. 
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Efficient Frontier forlBIVI-HWP 




Figure 17: Efficient frontier for a portfolio composed of two stocks: IBM and Hewlett-Packard. The dashed 
line represents the efficient frontier with respect to the second cumulant, i.e., the standard Markovitz efficient 
frontier, the dash-doted line represents the efficient frontier with respect to the fourth cumulant and the solid 
line the efficient frontier with respect to the sixth cumulant. The data set used covers the time interval from 
Jan. 1977 to Dec 2000. 
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Efficient Frontier forlBIVI-KO 
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Figure 18: Efficient frontier for a portfolio composed of two stocks: IBM and Coca-Cola. The dashed line 
represents the efficient frontier with respect to the second cumulant,i.e., the standard Markovitz efficient 
frontier, the dash-doted line represents the efficient frontier with respect to the fourth cumulant and the solid 
line the efficient frontier with repect to the sixth cumulant. The dataset used covers the time interval from 
Jan. 1970 to Dec 2000. 
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